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Application of: Teruhiko HAOIW ARA Confirmation No,; 1331 
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DOMAIN AVERAGING 

DECLARATION OF THE INVENTOR UND ER 37 C.RR, 8 1,131 

Assistant Commissioner for Patents 
Washington, DC 20231 

Sir 

I, Teruhiko HAOIWARA, do declare and state that: 

1 . I am a citizen of Japan residing at 941 3 Basson Drive, Houston, Texas, 77023. 

2. I am the sole named inventor of the invention disclosed and claimed in the above* 
identified patent application, Serial No. 09/803,819 filed on March 13, 200L 

3. I understand that U*S. Patent No. 6,452,389, having an effective date of February 7, 
2001, has been cited against the claims of the present application in an Office Action dated 
September 20, 2003. I am providing this Declaration to "swear behind" the effective date of 
the '389 patent. 

4. Attached hereto is Exhibit A, which is a copy of an invention disclosuie submitted by 
me to Halliburton Energy Services. The disclosure describes NMR logging using time 
domain averaging. 

5. I have reviewed the document in Exhibit A. I confirm that I conceived the invention 
disclosed in the document in Exhibit A before February 7, 2001 and diligently worked toward 
the filing of the present application on March 13, 2001. Although the dates have been 
redacted, I further confirm that the dates on the document are prior to February 7, 2001 . 
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6. I declare further diat all statements made in this Declaration of my knowledge are true 
and that all statements made on information and belief are believed to be true and further that 
these statements are made with the knowledge that willful false statements and the like so 
made are punishable by fine or imprisonment, or both, under Section 1001 of Title 18 of the 
United States Code and that such willful false statements may jeopardize the validity of the 
application or any patent issuing thereon. 



Exhibit A: Invention Disclosure • Time-domain average method to increase S/N ration in 
single event echo-train in high-speed NMR Logging and to sharpen the spatial 
resolution of NMR logging. 
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HALLIBURTON ENERGY SERVICES GROUP 

Page J_ of 4 INVENTION DISCLOSURE P.M. No. 

TITLE: Time-domain average method to increase S/N ratio in single event echo-train in high-speed NMR Logging and to 
sharpen the spatial resolution of NMR logging. 

FULL NAME OF SOLE or FIRST JOINT INVENTOR 

Teruhiko Hagiwara EMPLOYEE NO. G2966 ENTITY/DEPT.:JRAD 

First M. Last 

9415 Bassoon Drive Houston TX 77025 

Address Street City Stats zip 

FULL NAME OF SECOND JOINT INVENTOR 

EMPLOYEE NO ENTITY/DEPT.:_ 

Rrat M. Last 



Address: Street Qty State Zip 

FULL NAME OF THIRD JOINT INVENTOR 

EMPLOYEE NO. ENTITY/DEPT.: 

Ftfst M. Last 



Address; Street City State Zip 



Description of InvenUon: Give a brief description of the invention, and where possible a sketch. Do not write on back of this sheet If additional space is necessary, use additional disclosure sheets 
(R&S 231 1) and number consecutively. Sign and date each sheet and have each witnessed. Witnesses must be persons capable of understanding invention. If invented prior to date on this sheet 
give full details-date, to whom disclosed, identify drawings, if any, etc. If used In the field, specify for what company and designate vrall and date. 

This patent concerns a new processing method to increase signal-to-noise ratio (S/N) of NMR echo^ain logging data when NMR logging was 
perfomned at high logging speed and conventional stacking technique over multiple echo events cannot be applied to increase the S/N ratio. The 
new method enables us to run NMR logging at higher speed and obtain T2 inversion data without stacking echo train data from multiple events. 

Background 

One of the key information sought from the NMR logging data is the T2 decay time distribution which may be closely related to the pore size 
distribution of the rock formation. From the T2 decay time distribution, the capillary bound water volume may be estimated, for instance. The T2 
distribution is obtained by inversion of echo train data. Inversion of echo train is greatly affected by noise. In order to reduce the effect of noise 
and increase the signat-to-noise ratio, stacking of multiple echo train data Is commonly exercised, by assuming that noise is of random nature. 
For instance, Flg.1 illustrates that the average echo train from 10 events Is affected with much less noise than the echo train from single event. 

In actual logging, the tool Is moving at certain speed, and echo train from each event Is not collected at the same logging location or depth. When 
the logging speed is slow, the tool moves only short distance. Hence, one may assume that the earth formation Is of neariy constant property 
over the short distance, and stacking of multiple events from different depth may be acceptable. 

However, the logging speed is fast, the tool moves longer distance and the fomiation property may change significantly over the distance. As a 
result, stacking of echo train data from multiple events reduces apparent spatial resolution of measurements beyond the intrinsic resolution of the 
tool. 

Merit of Present Invention 

As we describe below, the present Invention uses time-domain average of single echo-train data to increase the S/N ratio, but not stacking of 
echo-train from multiple events, and hence, achieves higher spatial resolution. The S/N ratio can be increased also by application of both the 
time-domaln-average and stacking together. In this case, one needs less number of stacking events and higher resolution will be resulted. 

The concept of time-domain averaging is not new. For Instance, the US Patent No. 5,219,137 (Freedman) and 5,381,092 (Freedman) also utilize 
the time-domain averaging in order to compress logging data at downhole (In a wellbore) prior to transfer the data to a surface apparatus. 
However, the present invention uses the time-domain averaging, not for data compression, but to increase S/N and to enhance vertical 
resolution of the log measurements. In the present invention, ttie time-domain averaging is applied to tiie log data after data ti-ansmission. 

Present Invention 

The T2 decay time distiibution is obtained by inverting ttie time-domain echo ti^in data by assuming the physical model: 

Echo{t) = X^(^2 ) exp(-^ / ^2 ) + ^oise 
h 

where (p(T2) is tiie porosity (or population) of ttie pores corresponding to ttie exponential decay time T2. 

For simplicity, we omit ttie effect of partial polarization due to ttie longitudinal relaxation time Ti. It can be included, if necessary. One such 
example of echo-train is plotted in Fig.1 , where ttie echo train is generated from ttie T2 distiibution, (p{T2), of Flg.2. Because of noise, ttie inverted 
T2 distiibution obtained from ttiis single event echo-tialn. which is also shown in Flg.2, is different from ttie model input. 
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TITLE: Time-domain average method to increase S/N ratio In single event echo-train in high-speed NMR Logging and to 
sharpen the spatial resolution of NMR logging. 

To reduce the effect of random noise, it is commonly practiced to stack echo-trains from multiple number of events. See Fig.3. Fig.4 shows such 
stacked (averaged) echo trains obtained from 10 events. Significant reduction of random noise is observed. The inverted T2 distribution from the 
10 event stacked echo train is closer to the model input, as illustrated In Fjg.5. 

Stacking of multiple event data causes worsening of spatial resolution of NMR logging, as the tool is moving in usual logging measurements. 
The tool moves considerable distance at a high logging speed. The fonnation property may change significantly over the distance. While single 
event echo train provides infomnation of fomnation property for an interval determined by the intrinsic spatial resolution of the tool, the stacked 
echo train data provides formation property averaged over the distance to cover such multiple events. As a result, the spatial resolution of the 
measurement is greatly reduced beyond the intrinsic resolution of the tool. 

The S/N ratio can be increased without stacking echo train data from multiple events and thus without worsening the spatial resolution at high 
logging speeds, but by taking time-domain-average, (tn this method, the spatial resolution is not worsened at the price of the resolution in time 
domain.) See Fig.6. Namely, construct a time domain averaged echo train as follows: For time t, take an average of echo over the time interval 
of A, as defined by 

Echo^{t)= \dfEcho{f)l^. 
t 

It is crucial that the noise in the averaged echo is still random. To ensure that, the averaged echo separated by time interval of A is considered. 

Namely, the averaged echo train at t=to, to+A, to+2A, to+NA is used to obtain T2 distribution. Fig .7 shows such a time-domain averaged echo 

train averaged over 10 sampling points. In this example, the single echo train has 300 data points at trrE=1, 300 (with Te=1 .2 ms). As A=12 
ms (10 sampling points), the averaged echo has 30 (=300/10) data points. The inverted T2 distribution obtained from this time-domain-averaged 
echo is closer to the model input than the one from a single echo, as illustrated in Fig. 8. Note that the following equation is used to estimate J 2 
distribution: 

Echo^ (0 = Y,^{T^ F2 (1 - exp(-A / )) exp(-r / ) + Noise 
h 

= Z ) exp(-/ / ) + Noise 

h 

where 

Note that this T2 distribution mapped from this new echo train data has the resolution of A, in contrast to the resolution of unit time inten/al from 
tiie original echo train data. This is because the time domain averaging is perfomied uniformly over tiie entire time domain using the constant 
time interval A. This requirement can be relaxed in a more general time domain averaging method below. 

In order to obtain the T2 distribution in fiill range, especially at shorter T2, and to retain high T2 resolution at shorter T2, it is necessary to keep 
shorter averaging interval A at the eariier echo time tiiough the averaging interval has to be larger at later echo time where S/N ratio decreases 
considerably. Hence, tiie above window averaging can be modified for variable time interval Ak at time tk. Namely, we average tiie echo train at 
time tk over the window interval of Ak as, 

Echo^ {t, ) = \df Echo{f ) / A, 

It is again cmcial that the noise in the averaged echo is still random. To ensure ttiat, tiie averaged echo train sampled at t=to, ti= to+Ao+1 , t2= 
ti+Ai+1 , tN.i+AN-i+1 should be used to obtain T2 disti"ibution. Namely, tiie averaging interval should not overiap. The interval width Ak can be unit 
interval, consisting of single echo time. Fig.9 illusti^tes such time-domain averaging intervals, where the first five echo data are raw data without 
averaging, and tiie following 15 echo windows are two time unit lengtii for k=6-10, 4 unit length for k=1 1-15, 7 unit lengtii for k=16-20, and the 
rest is averaged over constant 10 unit length. The T2 distribution can be obtained from this time-domain-averaged echo by tine following equation: 

Echo,{t,) = X^(^2)7i(l - exp(-A, /r^WexpH, IT,) + Noise 
= Y.^(T,;A, (t, )) exp(-/, IT^)^ Noise 

where 

|(7;;A,(/,)) = ^(r,)r,(l-exp(-A,(^)/r,))/A,(/,) 
Note that tiie window widtti An depends on the echo time tk. 
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TITLE: Tlme-domaln average method to Increase S/N ratio in single event echo-train In high-speed NMR Logging and to 
sharpen the spatial resolution of NMR logging. 

Stacking of multiple event echo train data and time-domain averaging can be combined together. For instance, instead of taking 10- event 
stacking or time-domain averaging over 10 sampling points, one can take 5-event stacking and the apply time-domain averaging to increase the 
S/N ratio but not to worsen the vertical resolution much. 

Figures 




Figl. Single event echo-train 



Fig.2 Inverted Tj distribution from single echo-train data. 



One-event echo 



tntrinsic vartical retolutitin 
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Mutti-event echo 

apparent varticsl resolution 
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Rg.3 Multi-event echo stacking 





Rg.4 Started echo-train from 10 events 

One-event echo 



itrinsic vertical resolution 



ng.5 Inverted T; distribution from stacked echo-train data. 



Time-domain average 
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Rg,6 TimeKlomain average 
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TITLE: Time-domain average method to increase S/N ratio in single event echo-train in high-speed NMR Logging and to 
sharpen the spatial resolution of NIVIR logging. 




Fig.7 Time-domain averaged echtHrain Rg.8 Inverted T2 distribution from time-domain averaged echo-train 

Rg.9 Trme-domain averaged echchtrain with variable window vindtti 
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